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The azimuthal cos (j> and cos 2(j) modulations of the distribution of hadrons produced in unpolar- 
ized semi-inclusive deep-inelastic scattering of electrons and positrons off hydrogen and deuterium 
targets have been measured in the Hermes experiment. For the first time these modulations were 
determined in a four-dimensional kinematic space for positively and negatively charged pions and 
kaons separately, as well as for unidentified hadrons. These azimuthal dependences are sensitive 
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to the transverse motion and polarization of the quarks within the nucleon via, e.g., the Cahn, 
Boer-Mulders and Collins effects. 

PACS numbers: 13.88.-|-e, 13.60.-r 

Keywords: Semi-inclusive deep-inelastic scattering, azimuthal modulations, intrinsic transverse momentum 
and spin. 



I. INTRODUCTION 

Since the late 1960s the quark-parton model P, 0] has 
been used to describe the structure of the nucleon in 
terms of fundamental constituents. Their behavior inside 
the nucleon was parametrized in terms of parton distri- 
bution functions (PDFs). Historically PDFs depended 
only on the fractional quark momentum longitudinal to 
the nucleon direction of motion, x, and on the scale at 
which the distributions were probed, Q^, while trans- 
verse degrees of freedom were neglected. These PDFs 
have provided a good description of processes in which 
transverse spin and momentum are integrated over . 
However, transverse degrees of freedom are not a pri- 
ori negligible and are needed for a complete description 
of the nucleon. To account for transverse motion the 
PDFs have been generalized to transverse-momentum- 
dependent PDFs, known also as TMDs 

Already in the early days of the parton model it was 
realized that the inclusion of quark intrinsic transverse 
momentum, pr, leads to modifications of the cross sec- 
tions in high-energy reactions involving hadrons in the 
initial state, e.g . , in lepton-nueleon deep-inelastic scat- 
tering (DIS) [ij, IS- In particular, in semi- inclusive DIS 
transverse momenta give rise to azimuthal dependences 
of the distribution of the produced hadrons about the 
direction of the virtual photon 0, i, In 1978, 

R. Cahn discussed the emergence of cosine modulations 
in the semi-inelusive DIS eross section in the presence of 
non-vanishing transverse parton momentum using sim- 
ple kinematic considerations (Cahn effect) [UllS- Once 
included, interplay between the parton transverse mo- 
mentum and the partons's and nucleon's spins can gener- 
ate further azimuthal asymmetries, as, e.g., in the Boer- 
Mulders effect. The Boer-Mulders mechanism was in- 
troduced for the first time in 1997 [1] in relation to 
naive- T-odd effectfO |l9l . |20| . For a long time naive- T- 
odd effects were believed to vanish due to time-reversal 
invariance [2l|. Recently, it was shown that final and 
initial-state interactions can produce naive- T-odd effects 
without violating T- invariance [23 - |2^ . Because it in- 
volves only the parton spin and not the nucleon spin. 
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the Boer-Mulders mechanism is a good example of how 
spin-related effects may play an important role, even in 
unpolarized reactions. Measurements of these novel cor- 
relations provide insights into the so far poorly explored 
partonic transverse degrees of freedom and can be used 
to gather information, in a model-dependent way, about 
the elusive parton orbital motion. 

In DIS the structure of the nucleon is probed by the 
interaction of a high-energy lepton (Z) with a target nu- 
cleon [N] via the exchange of electroweak bosons. In 
the kinematic region accessed at Hermes it is a good 
approximation to consider onl y t he exchange of a single 
photon (Born approximation) |25| . In semi- inclusive DIS 
measurements, at least one of the hadrons {h) produced 
in the collision is detected in coincidence with the scat- 
tered lepton {V): 



I + N ^ I' + h + X, 



(1) 



where X represents the remaining, unobserved, final 
state. The polarization-averaged semi-inclusive DIS cross 
section can be written in a mo del- independent way by 
means of four structure functions 12a, ^3 : 



da,,,, = 



2-K 



dx dy dz dP,f j 



y 



xyQ^ 2(1 - e) 



2x 



{Fu 



+ y/2e{l + e) F™^ * cos + e F^'^'f' cos 20}. (2) 

Here, —Q^ is the squared four-momentum carried by the 
virtual photon. In the target rest frame, y is the frac- 
tion of the beam energy carried by the virtual photon 
and z is the fraction of the virtual photon energy carried 
by the produced hadron. The hadron momentum compo- 
nent transverse to the virtual photon direction is denoted 
Ph±, and (j) is the azimuthal angle of the hadron pro- 
duction plane around the virtual photon direction with 
respect to the lepton scattering plane (see Fig. [1]). The 
quantity a is the electromagnetic coupling constant, and 
7 — 2Mx/Q with M the proton mass. The structure 
functions Fuu,t, Fuu.l, F^T"^, F^°J'^'^ depend on x, Q^, 
z and Ph±', the subscript uu stands for unpolarized beam 
and target, while t (l) indicates transverse (longitudinal) 
polarization of the virtual photon, and e is the ratio of 
longitudinal to transverse photon flux. 

Besides intrinsic transverse parton momentum, 
perturbative-QCD (pQCD) effects, like gluon radia- 
tion, will also lead to azimuthal dependences in the 
semi-inclusive DIS cross section j29l [30l |. However, 
they contribute mainly at large values of Phi^, and are 
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FIG. 1. Depiction of the azimuthal angle 4> between the scat- 
tering plane, spanned by the three-momenta (Z, V) of incom- 
ing and outgoing leptons and the hadron plane, defined by 
the respective three-momenta of the virtual photon and the 
produced hadron, q and Ph, defined according to the Trento 
convention 28] . 



next-to-leading order in the strong coupling constant. 
For hadron transverse momenta that are small compared 
to the hard scale Q {Ph± <C Q), TMD factoriza- 
tion 0, H, allows for an expansion of the structure 
functions in powers oi 1/Q and expresses them in terms 
of convolutions of transverse-momentum dependent 
distribution and fragmentation functions. 

The TMDs parametrize the nucleon structure and the 
fragmentation functions describe how the struck parton 
evolves into the observed hadronic final state. For sim- 
plicity, in the following the weak dependence of the 
TMDs is not explicitly written. Each convolution can be 
classified according to the suppression, in powers of 1/Q, 
at which they contribute to the structure function. Not 
all contributions from all possible suppression levels have 
been calculated yet. In this work primarily contributions 
up to a suppression of (1/Q) will be considered, but, 
considering the low average attainable at Hermes, 
contributions suppressed as (1/Q)^ or higher may be 
not negligible. 

The structure function related to the cos 20 amplitude, 
pco^24>^ receives a single unsuppressed contribution: 



pcos20 



E 

9 



(3) 



Additional contributions arc present only at a suppres- 
sion of cx (l/Q)^, or higher. The sum symbol, stands 
for a quark-charge-squared weighted sum over quark fla- 
vors. The symbol (Kiyvi represents a weighte43 convolu- 
tion integral over the intrinsic momentum px and over 



2 The weights Wi (Eq. W2 (Eq. gj, and W3 and W4 (Eq. [Sj 
are kinematic factors depending on pT and fcy ; for their complete 
expressions see Ref. [27| . 



kx, the momentum transverse to the struck quark direc- 
tion that the hadron acquires during the fragmentation 
process. 

The contribution shown in Eq. [3] is called the Boer- 
Mulders-Collins effect (also often referred to as sim- 
ply the Boer-Mulders effect). It involves the Boer- 
Mulders distribution function hf''^{x,p'^) @, which de- 
scribes the correlation between the transverse polariza- 
tion and transverse momentum of quarks in an unpo- 
larized nucleon, and the Collins fragmentation function 
H^''^{z,kx) which describes the probability for a 

transversely polarized quark to fragment to an unpolar- 
ized hadron with a certain transverse-momentum direc- 
tion; both these functions are chiral-odd. As hard QED 
and hard QCD interactions preserve chirality, two chiral- 
odd functions need to appear in conjunction, in order to 
have a chiral-cven observable. 

As discussed above, naive- T-odd observables can be 
non-zero in conjunction with final- or initial-state inter- 
actions, which are reflected by the presence of non-trivial 
gauge links in the definition of the TMDs [H,!!!. This 
gauge link leads to the direct QCD prediction that naive- 
T-odd distribution functions must have opposite signs in 
semi- inclusive DIS and Drell-Yan reactions [2^. To date, 
this sign change has not yet been confirmed. 

There are no contributions to Fu°u'^'^ at a suppression 
l/Q- Not all contributions beyond a suppression of \/Q 
have been calculated, however a term 



)'E[/l(^'^'T)®W.i^?(^,4) 



CX 



(4) 



arising from the Cahn effect [l7|, ll8| is expected. The 
Cahn effect has recently received increasing attention, 
as it can provide information about the average trans- 
verse momentum of unpolarized quarks in unpolarized 
hadrons. It involves the convolution over transverse mo- 
menta of the spin-averaged distribution and fragmen- 
tation functions, fi{x,Px) and Dl{z,k'^), respectively. 
Their transverse- momentum- integrated correspondents, 
fi{x) and I?f (z), are well known [3ll, |32j. However, their 
Pt- and fcy- dependences are poorly constrained by mea- 
surements, and thus the convolution integrals in Eq. [3] 
can be estimated only approximately. Moreover, the av- 
erage intrinsic transverse momentum (px) may depend 
on the parton flavor; thus a flavor-dependent measure of 
the Cahn effect, via, e.g., semi-inclusive DIS of identified 
hadrons, is highly desirable. 

The first non-zero contributions to the structure func- 
tion F^"*^ in Eq. [2l which is related to a cos (f> ampli- 
tude, are suppressed as 1/Q; subsequent contributions 
are suppressed as (1/Q)'^. Among the various contribu- 
tions suppressed as 1/Q, several involve either a distri- 
bution or fragmentation function that relates to quark- 
gluon-quark correlations, and hence is interaction de- 
pendent and has no probabilistic interpretation. In the 
Wandzura-Wilczeck approximation 33| all these terms 
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are neglected, and only two contributions are considered: 
M 



pcosi, 
^ uu 



TABLE I. Numbers of charged hadrons for each data set (in 
millions). 



M 



(5) 



In the first line of Eq. [5] the Boer-Mulders-CoUins effect 
is recognizable, while, in the second line, the Cahn effect 
is present. 

Only a few measurements of cos 2(/) and cos (j) ampli- 
tudes in semi-inclusive DIS experiments have been pub- 
lished over the past 30 years |34| - [37| . Most measurements 
averaged over any possible flavor dependence as they re- 
fer to hadrons without type nor charge distinction, and 
only hydrogen target |34l 436| or hydrogen and deuterium 



targets combined together [37[ were available. Recently, 
the Clas collaboration measured non-zero cosine modu- 
lations for positive pions [ssj produced by semi-inclusive 
DIS off the proton. The Compass collaboration pre- 
sented preliminary cos 2(/) and cos 4> amplitudes in semi- 
inclusive DIS but has not yet published final results. 
In Drell-Yan experiments non-zero azimuthal modula- 
tions have been measured [40l - l45j that violate the Lam- 
Tung relation [i^. Such a violation can be ascribed to 
the Boer-Mulders distribution function, as pointed out 
in Ref. [13] . Sizable modulations have been extracted in 
pion-induced Drell-Yan reactions, where a valence quark 
and a valence antiquark annihilate. When a sea par- 
ton is involved, as in proton-induced Drell-Yan processes, 
the measured modulations become smaller, suggesting a 
small Boer-Mulders function for the sea. 

This paper presents cosine modulations for positively 
and negatively charged unidentified hadrons as well as 
for identified charged pions and kaons produced by DIS 
off hydrogen and deuterium targets. 



II. THE HERMES EXPERIMENT 

The cosine modulations described in the previous sec- 
tion were extracted from measurements performed at the 
fixed-target Hermes experiment. Hermes acquired data 
from 1995 to 2007 with various polarized and unpolar- 
ized gaseous targets internal to the Hera 27.6 GeV elec- 
tron/positron storage ring at Desy. In this paper results 
are presented that were extracted using only the pure 
hydrogen and deuterium targets, where the lepton beam 
scatters directly off neutrons and protons (with only neg- 
ligible nuclear effects in case of deuterium). The spec- 
trometer [4^ was a forward-angle instrument consisting 
of two symmetric halves above and below the horizontal 
plane defined by the lepton-beam pipe. Particles with 
polar angles within ±170 mrad in the horizontal direc- 
tion and between ±(40-140) mrad vertically could be de- 
tected. The collected data were processed with a tracking 
code involving event-level fitting based on a Kalman-filter 
algorithm [49[, which corrects the tracking parameters 
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for the effects from magnetic fields and accounts for all 
detector materials and known mis-alignments. 

Lepton-hadron separation with an efficiency better 
than 98% was achieved using the combination of several 
detectors: a transition-radiation detector [i^, a dual- 
radiator ring-imaging Cherenkov (RICH) detector [HO, 
[5l| . a lead and scintillator preshower detector, and a 
lead-glass calorimeter [52|. Hadron identification was 
performed using the RICH detector, taking into account 
the entire event topology simultaneously, rather than a 
single particle at a time. This provides improved par- 
ticle identification compared to earlier algorithms (see 
appendix IbI for further details). 



III. DATA ANALYSIS 



A. Data selection 



The data used for this work were collected during the 
2000-2007 periods, with both lepton beam charges and 
unpolarized hydrogen and deuterium targets, as summa- 
rized in table HI In order to guarantee high-quality data, 
each event had to meet several criteria, such as good per- 
formance of the particle identification and tracking detec- 
tors. Each selected track satisfies geometric constraints 
to ensure that it originated from the beam-target interac- 
tion region and also remained well within the acceptance 
of the spectrometer. 

Events with at least one lepton and one hadron de- 
tected in coincidence are included in the semi-inclusive 
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DIS event sample if they satisfy the foUowing kinematic 
requirements. The DIS region is defined here by the kine- 
matic constraints > 1 GeV^ and > 10 GeV^, 
where is the squared invariant mass of the initial 
system of virtual photon and target nucleon. As a conse- 
quence of these requirements and the limited angular ac- 
ceptance of Hermes, x and y are restricted to the ranges 
0.023 < a; < 0.6 and 0.2 < y. In addition the restriction 
y < 0.85 is applied, dictated by the energy threshold in 
the calorimeter to ensure high trigger efficiency. In or- 
der to suppress hadrons not originating from the struck 
quark (i.e., to suppress those from the target fragmen- 
tation region), the requirements z > 0.2 and xp > 0.2 
are applied. Here, xp = 2p^/\/s is the Feynman scaling 
variable, where is the hadron momentum component 
parallel to the virtual photon, and is the total energy 
in the 7*p center-of-mass system. 

To ensure good identification of hadrons by the 
RICH, identified pions are required to have momenta 
within 1 GeV < Ph < 15 GeV and kaons within 
2 GeV < Ph < 15 GeV. RICH weights are assigned to 
each hadron. These weights correspond to the probabil- 
ities that the hadron is a pion or a kaon, as determined 
from the RICH hadron type hypothesis and the identi- 
fication efficiency, computed from a Monte Carlo simu- 
lation of the RICH detector, which has been tuned to 
data. No RICH identification or RICH weights are ap- 
plied for the data sample of unidentified hadrons. To be 
consistent with the pion sample, the momentum restric- 
tion 1 GeV < P;i < 15 GeV is also applied to unidentified 
hadrons. In the calculation of all kinematic quantities 
that require particle masses, the pion mass was used for 
unidentified hadrons. This is as a good approximation as 
70% (88%) of positive (negative) hadrons are pions (see 
table U). 

The selected event sample is corrected for contamina- 
tion from leptons that do not originate from the scattered 
beam but rather come from lepton-pair production in de- 
tector material or meson Dalitz decay (Tr^/yy — ^ 7e~'"e~). 
These events amount to less than 1% of the total number 
of events, and are typically concentrated toward the high- 
y region, where their contribution reaches 8%. These con- 
taminating processes are charge symmetric. Therefore 
events passing DIS selection but with the wrong lepton 
charge constitute a control sample which is kinematically 
matched to the background events wrongly included in 
the semi-inclusive DIS sample. The correction is per- 
formed by assigning a negative weight to the events with 
lepton charge opposite to that of the beam. 



B. Extraction procedure 

Experimentally, the azimuthal modulations of the un- 
polarized cross section can be accessed via the cos n(j)- 



moments {n = 1,2) 



J^^ cos n0 dauu d(/) 



(6) 



where da^u is defined in Eq. [2l The moments are re- 
lated to the structure functions of interest via the 0- 
independent part of the cross section 



2 (cos ( 



v/2e(l + e) 
pcos20 _ 2(cos2^ 



(Fuu.T + eFuu.L), (7) 



(7(7, T ~ UU,L 



eF, 



(8) 



Extracting the cosine modulations of the unpolarized 
cross section from data requires disentangling them from 
a number of experimental sources of azimuthal modula- 
tions. At Hermes, due to the separation of the spec- 
trometer into two symmetric top-bottom halves, the az- 
imuthal acceptance is non-uniform. In addition, the ob- 
served kinematic conditions of each event may differ from 
the Born conditions at the hard electromagnetic vertex 
due to both physical and experimental effects. Events 
may be reconstructed with altered kinematic conditions 
due to external bremsstrahlung and multiple scattering 
in the detector material, or due to initial or final state 
radiation from the beam lepton (higher-order QED ef- 
fects). All of these effects lead to a miscalculation of the 
kinematic conditions, and can induce false cos n(j) modu- 
lations. 

To correct the data for kinematic smearing and false 
cosine modulations, a binned unfolding procedure was 
applied. A large Pythia6 [H^I Monte Carlo simulation 
(with approximately 20 times more events than the ex- 
perimental data) was generated, which uses the Jet- 
SET [HI] fragmentation model tuned to Hermes kine- 
matic conditions [5^. This simulation includes QED 
radiative effects calculated with Radgen [5^, and a 
complete GeantS [H^] simulation of the Hermes spec- 
trometer. All known instrumental and reconstruction ef- 
fects are simulated, including particle interactions with 
detector materials and detector responses that account 
for known inefficiencies. The simulated events then 
pass through the Hermes event-reconstruction algo- 
rithm, mimicking any possible tracking bias or ineffi- 
ciency present in the data. The simulated semi-inclusive 
DIS sample provides information on both the Born-level 
and the observed smeared kinematic conditions, and thus 
it can be used to build a matrix that describes the migra- 
tion of events between kinematic bins. This simulation 
is also used to define the events that smear into the mea- 
surement from outside the accepted kinematic range, and 
thus represent background events in each kinematic bin. 

A smearing matrix is constructed by normalizing the 
migration matrix to the Born cross section, taken from 
a Born-level PythiaG production, in each bin. As a re- 
sult of this normalization, the smearing matrix is a rela- 
tive quantity, and reflects the fraction of events that are 
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within the Hermes acceptance and their bin-by-bin mi- 
gration. In the limit of infinitely small bins, the smearing 
matrix is independent of the cross section model used to 
build it. However, the kinematic distribution of the back- 
ground events depends on the Born cross section model 
used to describe events generated outside the Hermes 
acceptance. The model dependence of the smearing ma- 
trix (due to finite sized bins) and the background are 
discussed in section [TVBl 

The simulated samples are normalized relative to the 
data via the inclusive DIS cross section, as determined 
by the LEPTO Monte Carlo generator To cor- 

rect the data for events smeared into the acceptance, the 
background events are subtracted from the normalized 
yields. The smearing matrix is then used to unfold the 
background-subtracted data, correcting for QED radia- 
tive effects, detector effects, acceptance and all the de- 
tector inefficiencies included in the simulation. 

The functional form 

^-f 6 cos 04- C cos 20 (9) 

is fit to the azimuthal distribution of the unfolded yields 
to extract the cosine modulations {cos (j))uu — B/2A and 
(cos20)„„ = C/2A. As the unfolding procedure is a lin- 
ear operation, it can be combined with the linear oper- 
ation of fitting. This was done with linear regression, 
where a was formed: 

^2 ^ (^data _ sX[3 fC-\a'^''^'' ~ SXp). (10) 

Here, cr'^'^*^ is the measured, background-subtracted 
yield, S is the smearing matrix, and C is a covariance ma- 
trix that includes the statistical uncertainties of data and 
background, and the statistical precision of the Monte 
Carlo used to construct the smearing matrix. The prod- 
uct XjS gives the fit function of Eq. [9] representing the 
Born-level event yield, with (3 the vector of parameters 
{A, B, and C) and X the block diagonal design matrix 
that includes a constant term (equal to 1) and the values 
of cos(0) and cos 2(0) for each bin. See appendix lA II 
for a detailed explanation of these matrices. 

The running conditions at Hermes changed from year 
to year. Thus the data of each year must be indepen- 
dently unfolded with the proper Monte Carlo production, 
including the appropriate experimental configuration. To 
combine the results from the various data sets, the for- 
malism of Eq. [To] has been extended to a procedure that 
at the same time unfolds each data set independently and 
fits the Born-level yields from all years simultaneously. 
Technical details about the full procedure are provided 
in appendix I A 11 

As the Born cross section depends on five kinematic 
variables (see Eq. [2|), this procedure is carried out on 
a five dimensional grid of kinematic bins. An analysis 
in fewer dimensions would implicitly integrate over vari- 
ables and mix together physics and experimental effects. 
This mixing can obscure the true signal, as demonstrated 
by Monte Carlo tests. A Monte Carlo simulation with 



TABLE II. Kinematic bin boundaries 



X 


0.023 


0.042 


0.078 


0.145 


0.27 


0.6 




y 


0.2 


0.3 


0.45 


0.6 


0.7 


0.85 




z 


0.2 


0.3 


0.4 


0.5 


0.6 


0.75 


1.0 


Ph± [GeV] 


0.05 


0.2 


0.35 


0.5 


0.7 


1.0 


1.3 





12 equi 


distant bins 











an isotropic distribution in at the Born level was run 
through the detector simulation. When unfolded in less 
than five dimensions, false modulations were extracted, 
which were of similar size as the physical moments seen 
in the data [59| . 

The binning used is reported in table HIl After the fit to 
the 0-dependence, the final four-dimensional (4Z)) cosine 
modulations represent fully differential results. Due to 
the unfolding procedure the results in the various kine- 
matic bins are statistically correlated as well as the re- 
sults for the (cos 0) and (cos 20) moments in each bin due 
to the fitting procedure. Therefore, the complete covari- 
ance matrix must be considered to avoid overestimating 
the statistical uncertainties in results projected on fewer 
dimensions. 



IV. SYSTEMATIC UNCERTAINTIES 

This section discusses the systematic uncertainties re- 
lated to the imperfect treatment of instrumental bias and 
inefficiencies by the smearing matrices. Systematic con- 
tributions related to residual model dependence of the 
unfolding procedure due to finite bin sizes are also eval- 
uated. 

In contrast to the case of unidentified hadrons, the sys- 
tematic uncertainties for identified pious and kaons in- 
clude an additional contribution from the RICH identifi- 
cation. In these cases, the Monte Carlo studies described 
below include a full simulation of the RICH detector. 



A. Instrumental effects 

The geometric acceptance of the Hermes spectrom- 
eter produces cosine modulations larger than the mea- 
sured signals; therefore a number of systematic checks 
have been performed and are listed in this section. 

The experimental apparatus experienced several ma- 
jor changes over time. The lepton beam charge changed 
as Hera alternated between accelerating electrons and 
positrons. For the last two years of data taking, the 
target cell was shifted in the beam-line direction closer 
to the forward spectrometer, and its length, initially of 
40 cm, was reduced by a factor of two. Different mag- 
netic fields were active in the target region in different 
years. The cosine modulations were extracted from data 
collected without any target magnetic field and with Ion- 
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gitudinal magnetic fields: a solenoid of 0.3 Tm strength, 
employed for a longitudinally polarized target in 2000, 
and the 1.0 Tm solenoid of the Hermes recoil detector 
installed after 2005. In addition, during shut-down peri- 
ods some detectors were moved in and out, and relative 
positions between detectors changed. All these altered 
conditions induce changes in the geometric acceptance; 
therefore each data taking period requires a dedicated 
simulation to properly correct for the acceptance. 

Despite those significant changes in running condi- 
tions, the cosine modulations extracted separately for 
each year are found to be mostly consistent. Small sys- 
tematic shifts between years are observed, which can be 
ascribed to effects not included in the simulations, and 
thus in the correction. These effects include residual de- 
tector misalignment not accounted for in the tracking 
algorithm, which are expected to change from one data- 
taking period to another. 

The time stability of the apparatus response was 
checked by measuring the azimuthal modulations gen- 
erated by the Hermes acceptance in short time intervals 
within the same data taking period. The test indicated 
that the azimuthal modulations of the acceptance are 
stable in time. 

The tiny instabilities (< 1% of the observed ampli- 
tudes) are highly dependent on the year under study, 
and, as the Monte Carlo does not simulate or correct for 
any of these instabilities, they can partially explain the 
small differences between cosine modulations extracted 
from different years. 

To take these differences into account, the signed dif- 
ference between moments extracted from each year and 
moments extracted from a simultaneous fit of the remain- 
ing periods was evaluated at the 4D-leve\. The modulus 
of the weighted average of these differences was added 
to the systematic uncertainty, and represents the largest 
contribution to the systematic uncertainty (^70% of the 
total systematic uncertainty). 

The extracted moments were checked for a sensitivity 
to a possible beam misalignment or slope with respect 
to the spectrometer axis, and misplacement of the spec- 
trometer dipole magnet. No significant effects have been 
found. The net beam polarization was found to be negli- 
gible. Additional instrumental sources that could gener- 
ate false azimuthal modulations have been tested by mea- 
suring sine modulations and cosine modulations higher 
than cos 20, which are not present in the unpolarized 
semi-inclusive DIS cross section in the single-photon ex- 
change approximation. No significant signals were found. 

The final moments discussed in section IVl and IVTl have 
not been corrected for possible binning effects or for 
RICH inefficiencies or cross-contaminations that were not 
accounted for (for identified pion and kaon samples only) , 
for example due to 0-dependence not accounted for in 
the RICH weights. The influence of these effects on the 
final moments was estimated by a Monte Carlo simula- 
tion. For each particle type under study, a AD model 
of the measured cosine modulations was extracted from 



the fully differential final moments by means of a AD 
parameterization (details in appendix lA 2|1 . Through 
an accept/reject procedure, those models were used to 
alter, at the Born-level, the underlying distribution in 
an originally azimuthally uniform PythiaG production 
that includes the full spectrometer simulation and QED 
radiative effects. For the identified hadron cases, both 
pion and kaon models were implemented, to account for 
cross-contaminations between modulations. The protons 
constitute the remaining significant part of the hadron 
sample, but a model was not extracted for them. The 
sensitivity of the test to proton modulations was checked 
by implementing a model for protons that was either 0- 
independent or with modulations identical to that of pi- 
ous. The proton model input to the Monte Carlo is found 
to have very little impact on the test results. 

This simulation, modeled to reproduce the measured 
cosine azimuthal modulations, was used as a surrogate for 
the data in the entire analysis procedure, and cosine mod- 
ulations for pions and kaons and unidentified hadrons 
were extracted. The extracted moments were found to 
agree with the input models, verifying that the unfolding 
algorithm is able to extract the implemented modulations 
after correcting for all instrumental and QED radiative 
effects included in the simulation. The small discrepan- 
cies between the extracted moments and the input model 
provide an estimate of systematic uncertainty due to the 
unfolding procedure, binning effects and RICH weights 
(in the case of the identified pion and kaon samples). 



B. Model dependence 

The unfolding procedure described in section fill Bl can 
be affected by two different sources of model dependence. 
The unfolding correction is based mainly on two objects: 
a smearing matrix, describing the migration of events 
between bins, and a background estimation, describing 
the events that are smeared into the kinematic bins from 
outside the geometric/kinematic acceptance. 

In a fully differential analysis and in the limit of in- 
finitely narrow bins, the smearing matrix is independent 
of the models underlying the Monte Carlo event genera- 
tor used to produce it. Residual model dependence due to 
finite bin sizes was tested by comparing data azimuthal 
moments extracted using smearing matrices computed 
with different models for the azimuthal dependent part of 
the cross section: the standard, (/)-independent, Pythia6 
cross section, and the altered PythiaB cross section that 
includes the AD cosine model extracted from data, as de- 
scribed in section HV Al As expected, no significant dif- 
ferences were observed in the extracted moments. 

To test the model dependence of the background a 
similar procedure was used, and cosine modulations ex- 
tracted with different models for background evaluation 
were compared. Again the two models used were the 
standard, 0-independent PythiaG cross section, and the 
PythiaB cross section modified to include the AD cosine 
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model, which was extrapolated into the unmeasured 
region not covered by the detector acceptance. The 
moments from data were found to be weakly sensitive 
to the azimuthal dependence of the background. The 
differences between the moments extracted with the two 
models were used to estimate the systematic uncertain- 
ties from the model dependence and were combined with 
the other systematic uncertainties in quadrature. 



C. Calculation of systematic uncertainties 



TABLE III. Kinematic ranges of integration. 



0.023 - 0.27 



0.042 - 0.27 



Kinematic Range A 
y z 
0.3 - 0.85 0.2 - 0.75 

Kinematic Range B 

y z 

0.3 - 0.7 0.2 - 0.6 



Ph± [GeV] 
0.05 - 1.0 



Ph± [GeV] 
0.2 - 0.7 



The systematic uncertainties described above are sub- 
ject to statistical fluctuations due to the finite statistical 
precision of the Monte Carlo simulations used to calculate 
them. To average out these statistical fluctuations, each 
systematic contribution was smoothed by fitting it to a 
AD linear function. Higher order 41? polynomials were 
tested, and provided final systematic uncertainties of sim- 
ilar size. The final systematic uncertainty was calculated 
by adding each smoothed contribution in quadrature. 



V. FULLY DIFFERENTIAL RESULTS 

The AD analysis described in section IIII Bl provides 
access to the full kinematic dependences. The final 
moments, in four dimensions, for positive and negative 
unidentified hadrons, pions, and kaons produced from hy- 
drogen and deuterium targets, are available online [GOj. 
These fully differential moments represent the complete 
set of results of this analysis and can be used to test 
theoretical models. 

The moments are accompanied by the covariance ma- 
trix describing the statistical correlations, the total sys- 
tematic uncertainties for each bin (all contributions are 
added in quadrature), and the average (x), (y), (z), 
{Ph±), and (Q^) values for each bin. It is not possible 
to make a measurement in every (x, y, z, Ph±) bin due 
to (a) kinematic constraints that exclude some portions 
of the AD space, (b) the not-uniform distribution of the 
underlying cross section across the rectangular kinematic 
binning, (c) the detector acceptance, and (d) the limited 
statistical precision of the data. Bins that do not con- 
tain enough events to make a measurement (and typically 
also have statistical uncertainties larger than unity) are 
denoted in the database with all moments and average 
kinematics set to zero. In the covariance matrix, ele- 
ments corresponding to such bins have diagonal element 
values of one and non-diagonal element values of zero. 

A visual representation of the bins where a measure- 
ment is not possible can be found online [6l| for the sta- 
tistically poorest data set for each particle type (pions, 
kaons. and all hadrons). This tool also allows the user to 
integrate the moments in an arbitrary kinematic range 
(following the procedure described in the next section). 



VI. RESULTS FOR FIXED KINEMATIC 
RANGES 

The fully differential moments provide the maximum 
information from this measurement. In order to gain 
a qualitative picture of the behavior of the moments, 
a projection to one dimension (ID) was performed by 
a weighted integration of the moments over three vari- 
ables, highlighting the dependence of the moments on 
the fourth variable. In order to achieve an integral over 
the selected kinematic ranges, the moment in each kine- 
matic bin is folded with the (/)-integrated unpolarized 
semi-inclusive DIS cross section in this bin, normalized 
to the same value integrated over the whole kinematic 
range of the projection. The necessary input cross section 
for this integration was extracted directly from Hermes 
data [ell. 

As anticipated in section |Vl it is not possible to make 
a measurement in every bin. Therefore, restricted kine- 
matic ranges of integration were chosen to minimize the 
number of bins in the integration where a measurement 
is not possible. In addition, the z-bin from 0.75 to 1.0 
is excluded as a large fraction of events in this kinemat- 
ics range contain decay products of exclusively produced 
hadrons, for which standard factorization might be bro- 
ken. Hadron and pion results are integrated over the 
kinematic range A, given in table IIIll This table also lists 
the reduced kinematic range B used for kaons, which have 
comparatively lower statistical precision. The average 
kinematic values for each integrated bin for the ranges 
A and B are shown in Figs. [5] and [31 respectively. The 
strong correlation between x and (and thus between 
X and y) shown in the bottom left subplot is due to the 
Hermes acceptance. The other kinematic variables are 
weakly correlated. 

In the integration, bins with large statistical uncertain- 
ties (larger than unity) are not included as they do not 
provide meaningful information to the integral and inflate 
its uncertainty. The effect of excluding these bins is esti- 
mated and included in the systematic uncertainty of the 
projected results, as described at the end of this section. 
In figures that compare results from various hadron types 
or charges, the ranges of integration have been chosen so 
that only the bins that provide a measurement in the 
statistically poorest data set are included. This assures 
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FIG. 2. Average kinematics for the integration range A of table IIIII as extracted from a An Monte Carlo (shown here for 
positive hadrons on hydrogen; other cases exhibit only minor deviations). 
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FIG. 3. As in Fig. [2] but for the kinematic range B of table Hill 
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that the data sets have identical integration regions and 
thus allows for the results to be compared in a consistent 
way. 

The integrated azimuthal modulations were found to 
be weakly sensitive to the semi-inclusive DIS cross sec- 
tion used for the integration. This sensitivity was as- 
sessed by using the cross section extracted from Hermes 
data as well as the cross sections implemented in two 
Monte Carlo simulations tuned to reproduce the Her- 
mes measured yields: Pythia6 [HI] and GMC_TRANS. 
The Hermes GMC_trans generator uses the Cteq6 dis- 
tribution functions and the DSS fragmentation func- 
tions [12,1111 ■^ith the P?i_L-dependence based on a Gaus- 
sian ansatz. In particular, the transverse momenta px 
have a non-constant z -dep endence as observed from a 
fit to Hermes data [il], while for Pythia6 this z- 
dependence is flat. The sensitivity to semi-inclusive DIS 
cross section used for the integration has been added to 
the systematic uncertainty. 

The systematic uncertainties for the results projected 
in ID are composed of the uncertainties discussed in 
section llVi the sensitivity to the cross section used in 
the integration described in the previous paragraph, plus 
the additional uncertainty added by excluding some bins 
from the projection. The Monte Carlo production mod- 
ified to reproduce the measured azimuthal distribution 
(described in section lTV A[) was used to evaluate the effect 
of the bins excluded from the integration. The difference 
between including or excluding these bins in the integra- 
tion of the simulated moments was added to the other 
systematic contributions. Each systematic contribution 
was independently projected onto the single variable be- 
fore the smoothing described in section PV CI After the 
projection, the systematic contributions were smoothed 
with a ID linear fit and then added in quadrature. 



A. Results for charged pions 

The cosine modulations for charged pions, projected in 
the kinematic range A (table IIIip , are presented in this 
section. All pion samples are projected only including 
bins that provide a measurement in every data sample, 
which restricts the integration to those bins with a mea- 
surement in the statistically poorest data sample, i.e., the 
sample for negative pions produced from the deuterium 
target. 



1. Pion cos2(^ amplitudes 

Figure m shows the cos 20 amplitudes 2{cos2(t))uu for 
pions extracted from hydrogen and deuterium data, pro- 
jected versus x, y, z, and Phi.- Different magnitudes and 
opposite signs of the amplitudes are observed for oppo- 
sitely charged pions. In particular, positive cos 20 ampli- 
tudes are extracted for negative pions, while for positive 
pions the moments are compatible with zero, but tend to 



be negative in some kinematic regions. The amplitudes 
for positive and negative pions also exhibit different kine- 
matic dependences. This is particularly evident in their 
dependence on z: in the integrated kinematic region pre- 
sented here, the magnitudes for positive pions have no 
clear kinematic dependence, while they rise with z for 
negative pions. The amplitudes increase in magnitude 
with P/ix for both 7r+ and 7r~, but with opposite signs. 

Up to a kinematic suppression of ( 1/Q) , the cos 20 am- 
plitudes only contains a single, unsuppressed term, the 
Boer-Mulders-Collins effect, i.e., the convolution of the 
Boer-Mulders distribution function hj-{x,p^) and the 
Collins fragmentation function H^{z,k^) discussed in 
section m For a hydrogen target, scattering off up quarks 
is expected to dominate the reaction (u-dominance) , both 
because the proton consists of more up quarks than 
down quarks, and because the elementary lepton-quark 
cross section is proportional to the squared quark electric 
charge (e^), which gives an additional factor of 4 for up 
quarks compared to down quarks. The Collins function 
was recently found to have a similar magnitude but op- 
posite sign for fragmentation of up quarks into positive 
{favored fragmentation) and negative pions {disfavored 
fragmentation) [65|-|69|. This would result in different 
signs for pions of opposite charge, which is in agreement 
with the data. 

The similarity between hydrogen and deuterium re- 
sults seems to indicate that the Boer-Mulders distribu- 
tion function has the same sign for up and down quarks, 
as shown in Ref. [T^I and in Ref. [7l|, and anticipated 
in Refs. [zl, Although they are similar, for posi- 

tive pions the deuterium results seem to be systemati- 
cally closer to zero with respect to the hydrogen results; 
this might be due to a different magnitude of the Boer- 
Mulders function for up and down quarks. 

Model calculations [t^ - Itt} of the contribution of the 
Boer-Mulders-Collins effect to the (cos 20) moment 
are in qualitative agreement with the moments reported 
here. In particular, the opposite sign for oppositely 
charged pions seems to be a signature of the Collins ef- 
fect. 

Equation [3] only includes terms up to a suppression 
(l/Q), but at a suppression of (1/Q)^ there is at least 
one additional term that includes the Cahn effect (see 
Eq. HI). The restricted range of the Hermes data 
does not allow for a conclusive study that disentangles 
the leading term from the suppressed terms. Nonethe- 
less, an attempt to describe preliminary Hermes results 
in a more complete way has been done in Ref. [tgI ] , where 
the authors evaluated this suppressed Cahn contribution 
to the cos 20 amplitude, assuming a flavor-blind Cahn 
term, i.e., a flavor-independent (p^). The comparison 
of this calculation to data indicates that, in the Her- 
mes kinematic regime, the Cahn term is smaller than ex- 
pected or is counteracted by additional terms that have 
been neglected. In the same paper, a possible Cahn 
flavor-dependence was also estimated by varying the (p^) 
for down quarks while maintaining a fixed (p'^) for up 
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quarks; no significant changes were observed in the cal- 
culated Cahn term. However, this test was performed on 
a hydrogen target, and not a deuterium target where the 
results might be more sensitive to the down quarks. 
In Ref. 
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the authors attempted to simultane- 
ously describe preliminary unidentified hadron cos 2^ 
amplitude extracted at Hermes and Compass [H]. 
The Boer-Mulders-Collins effect is described using the 
Collins fragmentation function from Ref. [66j while for 
the Boer-Mulders function, the same functional form 
that was used for the Sivers function [t^ was applied. 
In the calculation the Cahn effect is also included, 
which is sensitive to the quark average transverse mo- 
menta. The previously reported average momentum of 
0.25 GeV^ [iOl describes the Compass data well. 
In contrast, the Hermes data is better described by the 
lower value of (p^) ~ 0.18 GeV^, leading to a smaller 
Cahn effect at Hermes. This is in accordance to the 
broadening of the pt distribution when considering 
evolution, as observed in Ref. ISli ■ 



qualitatively describe the measured amplitudes for cos 0. 
To date, only one prediction has been published for the 
Hermes cos0 amplitude |82|, which includes the Cahn 
term only. The amplitudes predicted are larger than 
the measurements, suggesting that the Cahn contribu- 
tion at Hermes is smaller than expected. As in the 
case of (cos 20)^^, at least part of the discrepancy can 
be related to a (p^) that, in Hermes kinematic condi- 
tions, is smaller than 0.25 GeV^. Moreover, the modeled 
Cahn term cannot describe the observed difference be- 
tween TT"*" and 7r~, as it was assumed to be flavor-blind. 
This implies that for a qualitative description of the mea- 
sured cos4> amplitudes, a more complex Cahn contribu- 
tion, or additional flavor dependent contributions, like, 
e.g., the Bocr-Mulders-Collins effects, must be taken 
into account. Furthermore, in addition to the Cahn and 
the Boer-Mulders terms, the structure function F^}'^'^ 
includes four terms related to quark-giuon-quark corre- 
lators which have not been taken into account in this 
interpretation of the data, as little is known about the 
underlying physics. 



2. Pion cos (f} amplitudes 

The cos (/) amplitudes come suppressed as 1/Q in the 
hadron cross section, and, in contrast to the cos 20 ampli- 
tudes, several terms contribute at same level of suppres- 
sion (Eq. [5|). Results for the cosc/) amplitudes 2(cos0)t7t7 
extracted for pions from hydrogen and deuterium data 
are shown in Fig. [5] Results extracted from hydrogen 
and deuterium are similar, but deuterium results for pos- 
itive pions are smaller than hydrogen results. This could 
be related to flavor dependence of the contributions in- 
volved in the amplitudes. The cos (j) amplitudes are found 
to be negative for both positively and negatively charged 
pions, but for positive pions they are in general larger 
in magnitude. For both positive and negative pions, the 
magnitudes increase with the pion energy fraction z. 

The z dependence of the amplitudes can be interpreted 
in terms of the Cahn effect. Indeed, Cahn anticipated a 
rise of amplitudes with z due to the reduced dilution by 
the random transverse momentum that the pions acquire 
during fragmentation [13, [13 ■ At high z the amplitudes 
for oppositely charged pions are very similar and reach 
their largest magnitude (up to —0.2). Different behav- 
iors are observed for oppositely charged pions versus Ph±_ . 
The magnitude of the amplitudes for positive pions in- 
creases with Ph±, supporting the Cahn expectations of a 
signal proportional to transverse momentum. But, this 
trend is not observed for negative pions, suggesting that 
contributions to Eq. [5] other than Cahn possibly coun- 
terbalance the increase with Ph± . The Cahn term is ex- 
pected to be weakly sensitive to flavor, as discussed in 
the previous section. As in the case of the cos 20 am- 
plitudes, the difference between oppositely charged pions 
can be generated by flavor dependent contributions, like, 
e.g., the Boer-Mulders-Collins effect. 

In contrast to the cos 20 amplitudes, no model can 



B. Results for charged kaons 

This section presents the cosine modulations extracted 
for charged kaons, projected in the reduced kinematic 
range B of table IIIII All kaon samples are projected us- 
ing bins that provide a measurement in the negative kaon 
sample produced from the deuterium target, which is the 
statistically poorest kaon data sample. No model calcu- 
lation for kaons is available to date. 



1. Kaon cos 20 amplitudes 

The cos 20 amplitudes extracted for charged kaons are 
shown in Fig. [6] for the hydrogen and deuterium targets. 
The amplitudes are large in magnitude (up to —0.1), and 
have the same negative sign for both positive and nega- 
tive kaons, in contrast to the trends observed for pions. 
This may be interpreted by considering the kaon's quark 
content: the valence quark content of mesons is us, 
and therefore K'^ production is expected to receive a 
large contribution from lepton scattering off up quarks 
(u-dominance). A favored Collins fragmentation func- 
tion is expected to be involved in this case, as in the case 
of 7r+. In the framework of the Artru model [s^, all 
favored Collins functions describing fragmentation into 
spin-zero mesons have the same sign. Therefore, the 
Boer-Mulders-Collins effect for positive pions and kaons 
is expected to have the same sign, as observed in mea- 
surements. Nothing is known about the Collins frag- 
mentation function into kaons. A significant contribu- 
tion from sea quarks cannot be excluded. For example, 
strange quarks may contribute at a; < 0.1 [s^ . as sug- 
gested by unpolarized fragmentation, where the strange 
quark fragmentation function into K'^ appears to be sig- 



12 



3 
3 



0.1 



-e- 

CNJ 

<g 0.05 

o 

CM n 



-0.05 
Z> 0.1 

<5 0.05 

o 



-0.05 



71" 




■ e p ^ e 7t X 

° e d ^ e 7t X 



71 



i 



4 




10- 



0.4 



0.6 



0.8 

y 



0.4 



0.6 



0.2 0.4 



0.6 



0.8 1 

[GeV] 



FIG. 4. cos 2(j> amplitudes for positive (upper panels) and negative (lower panels) pions integrated over the kinematic range 
A of table [Till Closed and open squares are for amplitudes extracted from hydrogen and deuterium targets, respectively. The 
inner bar represents the statistical uncertainty; the outer bar is the total uncertainty, evaluated as the sum in quadrature of 
statistical and systematic uncertainties. Points have been slightly shifted horizontally for visibility. 



3 

0.1 




-0.1 
-0.2 



3 

0.1 




-0.1 
-0.2 



- 71+ : 






■ e p ^ e 7t X 

° e d ^ e 71 X 


'- ^ ^ i I ■ 


1 ... 1 ... 1 .■ 


9 , - 



* f ■ 

... 1 ... 1 ... 1 ... 1 ... r 


■ ^" f ■ 

1 


J 

* 1 f 

. 1 , , . 1 . , , 1 . 


--^-*-f 

- 1 . 

I- 


, , , 1 , . . 1 . . . 1 , , , 1 , , , 1' 



10-^ 



0.4 0.6 



0.8 

y 



0.4 0.6 



X y Z 

FIG. 5. As in Fig. |4l but for the cos (j> amplitudes. 



0.2 0.4 0.6 0.8 1 
> 

h_L 



P., [GeV] 



13 



nificantly larger than the fragmentation function for up 
quarks into [H]. A substantial difference between 
the strange and the up and down Collins fragmentation 
functions would play an important role in the observed 
moments. 

For negative kaons the situation is even more compli- 
cated, as its valence quark content (su) does not include 
any quarks in common with the valence structure of the 
target. Therefore, even larger contributions can be ex- 
pected to originate from the sea and from disfavored up 
quark fragmentation. 

Similar kaon amplitudes are extracted from hydrogen 
and deuterium targets. This may reflect similar contri- 
butions from u and d quarks, as well as a potentially 
substantial contribution from strange quark fragmenta- 
tion, which is expected to be the same for neutron and 
proton targets. Contrary to pions, the positive kaons 
show hydrogen results closer to zero, which might reflect 
a different magnitude for the Boer-Mulders function of 
different quark types, or the increased role of disfavored 
up quark fragmentation for proton targets. 

2. Kaon cos <f) amplitudes 

The cos (f> amplitudes for kaons are shown in Fig. [7| 
for hydrogen and deuterium targets. Large negative (up 
to —0.2) amplitudes are extracted for positive kaons, 
slightly rising with z and Ph±. The amplitudes are even 
larger in magnitude than those for positive pions, which 
suggests a large contribution from the Boer-Mulders- 
Collins effect, which was found to be large for K'^ in 
the previous section. Negative kaons instead show re- 
sults compatible with zero. The similarity between the 
cos 2(f) amplitudes for positive and negative kaons may 
mean that the Boer-Mulders-Collins effect is relatively 
insensitive to kaon charge. Thus, the significant differ- 
ence in the cos 4> amplitudes for positive and negative 
kaons points to either a flavor dependence of the Cahn 
contribution (e.g. from strange quarks) or a significant 
contribution from the interaction dependent terms that 
have been otherwise neglected in this discussion. Simi- 
lar results are extracted for scattering off hydrogen and 
deuterium. 



C. Results for unidentified charged hadrons 

In this section the cosine modulations extracted for 
unidentified hadrons and projected in kinematic range 
A (table are presented. As for identified charged 
hadrons, individual kinematic bins are included in the in- 
tegration only if they provide a measurement in the sta- 
tistically poorest unidentified hadron data sample, i.e., 
negative hadrons produced from a deuterium target. 

As the majority of the unidentified hadrons consists of 
pions (>70 — 88%, depending on the hadron's charge), 
the amplitudes of unidentified hadrons are very similar 



to those of pions, and most of the arguments from the 
discussion of the pion results also apply here. However, 
as no hadron identification was required, the systematic 
uncertainty for the unidentified hadron sample does not 
include a contribution from the RICH identification. The 
remaining hadrons are in large part kaons (^10%), and 
protons (^10%). As no theoretical model has evaluated 
the cosine modulations for kaons and protons, no predic- 
tions exist for the unidentified hadron sample. 



1. Hadron cos 2(/) amplitudes 

Figure [8] shows the cos 2(j) amplitudes of unidentified 
hadrons extracted from hydrogen and deuterium data. 
Different amplitudes are extracted for oppositely charged 
hadrons; in particular, they are of opposite sign, as in 
the case of pions. Similar amplitudes are observed for 
hadrons produced from hydrogen and deuterium targets. 



2. Hadron cos </> amplitudes 

Results for the cos amplitudes extracted from hy- 
drogen and deuterium data are shown in Fig. [9l They 
arc found to be negative for both positively and nega- 
tively charged hadrons, but larger in magnitude for the 
positive hadrons. Hadrons produced using hydrogen and 
deuterium targets result in similar amplitudes, but small 
differences can be observed for that refiect the behav- 
ior of the TT"*" amplitudes. 



D. Comparison of amplitudes for various hadron 
types 

In order to compare the cosine modulations extracted 
for the various hadron types, all samples were projected 
in the smaller integration range B of table IIIII and bins 
were only included if they provided a measurement in the 
negative kaon sample produced from a deuterium target, 
which is the statistically poorest data sample. FiguresfTOl 
\TT\ [T2l and [13] show the comparisons of the cos 2(j) and 
cos (f) amplitudes for the various hadron types produced 
on a hydrogen and deuterium target. The kaon moments 
are substantially larger in magnitude than those of the pi- 
ons, with the exception that the K~ (cos0)„„ moments 
are compatible with those of 7r~, although their large 
uncertainties also make them compatible with zero. The 
{cos2(f>)uu moments for negative kaons not only have a 
larger magnitude but also the opposite sign as the pion 
moments. A magnitude of amplitudes larger than 
that for 7r+ was already observed in the case of the am- 
plitudes measured in transverse-target single-spin asym- 
metries where the Collins fragmentation function cou- 
ples to the transversity distribution function (67j . The 
large amplitudes for kaons suggest a Collins effect that 
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FIG. 6. As in Fig. |4l but for charged kaon amplitudes integrated over the kinematic range B of table IIIll 
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is larger for kaons than for pions; in addition, the differ- 
ences with respect to pions can arise from a significant 
role of strange quarks in kaon production. The modula- 
tions extracted for unidentified hadrons and pions have 
similar trends, although some differences arc observed, 
particularly for the cos 2(j) amplitudes. The discrepancies 



between hadrons and pions are generally consistent with 
the observed kaon moments. 
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As in Fig. m but for unidentified charged hadrons. 
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E. Charge difference 

The systematic uncertainties of results in Figs. [3HS] are 
highly correlated for positive and negative hadrons of the 
same type, as they were measured under the same ex- 



perimental conditions. It is therefore useful to provide 
the difference between the amplitudes of positive and 
negative hadrons, where many systematic uncertainties 
cancel. The charge difference provides more strict con- 
straints for models, as it accounts for correlated system- 
atics between hadrons of the same type, but different 
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FIG. 12. As in Fig. [TOl but for a deuterium target. 



3 

(0 

o 

CM 



0.2 

0.1 



-0.1 

-0.2 

-0.3 
0.2 

0.1 



-0.1 

-0.2 

-0.3 











: A e d ^ e h X 
□ e d ^ e 7t X 
: oed ^eKX : 


'■■ p B ^ 

' ■ ■ 


- ^ B ^- 

: ^ 1 
I.I.I. 


1 ■ 


; i ] 
1 , 1 . 1 


: § ^ : 





3 

(0 

o 
CVJ 



4 



.4. 



I 



' ' ' ' I ' ' I ' I ' ' I ' ' I ' I ' ' I ' I ' 



10-1 0.4 0.5 0.6 0.7 

X y 



0.4 



0.5 



0.6 0.3 0.4 0.5 0.6 



FIG. 13. As in Fig. 1111 but for a deuterium target. 



charge. In addition some hadron-flavor blind contribu- 
tions to the moment may be suppressed, e.g., a Cahn 
effect as considered so far in most phenomenological ap- 
proaches. In that case, both the cos (p and cos 2(/) charge 
difference amphtudes are expected to have an increased 
sensitivity to the Boer-Mulders-Cohins effect. 



For each hadron type the charge difference of the re- 
spective amphtudes was evaluated, and its uncertainty 
was computed, taking into account the correlations. The 
results are shown in Figs. [TH [151 and [TBI for pions, kaons, 
and unidentified hadrons, respectively. For pions and 
unidentified hadrons the charge difference is significantly 
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FIG. 14. DifTerence between the amplitudes of negative and positive pions within the kinematic range A of table IIIII 
2(cos2(/>)^- — 2(cos2(^)^+ (upper panels), 2{cos (j})^- — 2(cos(^)^+ (lower panels). Closed and open symbols are for ampli- 
tudes extracted from hydrogen and deuterium targets, respectively. Points have been slightly shifted horizontally for visibility. 
Uncertainties as in Fig. 3] 
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non-zero over nearly the entire kinematic range. For 
kaons. a large charge difference is observed for the cos0 
amplitudes, while the difference for the cos 20 amplitudes 
is compatible with zero. The different behavior of kaons 
with respect to pions suggests an important contribution 
to the modulations from scattering off strange quarks, 
or, more generally, from scattering off the sea. 



VII. CONCLUSIONS 

Hermes measured fully-differential (4Z?) cosine modu- 
lations for charged pions, kaons, and unidentified hadrons 
produced in semi-inclusive DIS off unpolarized hydrogen 
and deuterium targets. In the TMD framework, these 
amplitudes can be interpreted by a non-zero intrinsic 
transverse momentum of quarks (Cahn effect) and by cor- 
relations between the quark's transverse polarization and 
its transverse momentum and the transverse momentum 
of the produced hadron (Boer-Mulders-Collins effect). 
However, considering the low average attainable at 
Hermes, contributions suppressed as (1/Q)^ and higher 
may be not negligible. 

To date, the cosine modulations presented here rep- 
resent the most complete data set on the subject, and 
allow access to flavor-dependent information on the in- 
ternal degrees of freedom of the nucleon. The extracted 
AD amplitudes, with their full covariance matrix, pro- 
vide the maximum information from this measurement, 
and can be used to guide model construction in a fully 
differential way. 



In addition, the amplitudes have been presented as 
projections over specific integration ranges of the four 
kinematic variables x, y, z, and Ph±- Significant differ- 
ences are observed for moments extracted for oppositely 
charged pions, interpreted as being due to the convolu- 
tion of the Boer-Mulders and the Collins functions. The 
extracted amplitudes for pions and kaons show differ- 
ent sizes and kinematic dependences. In particular, the 
cos 20 amplitudes for kaons are larger than for pions, 
and do not change sign for oppositely charged kaons. 
The former may be due to different features of Collins 
fragmentation into kaons and pions, possibly related to 
a significant contribution from strange quarks to kaon 
production. 

The amplitudes extracted from hydrogen and deu- 
terium targets are found to be similar, but slightly dif- 
ferent for positive pions and kaons. For pions, a similar 
size can be due to the Boer-Mulders distribution func- 
tions with the same sign for up and down quarks. For 
kaons this can be due to a similar contribution from up 
and down quarks along with a similar strange sea dis- 
tribution in protons and neutrons. The slight differences 
for hydrogen and deuterium targets for positive pions and 
kaons might be related to a slightly different magnitude of 
the Boer-Mulders functions for the different quark types. 
The difference of moments between positively and nega- 
tively charged hadrons is not compatible with zero for all 
hadron types except for the cos 20 kaons. The different 
behavior of kaons with respect to pions suggests a sig- 
nificant contribution to the modulations from scattering 
off strange quarks, or, more generally, from scattering off 
the sea, or from fragmentation of light quarks into kaons. 
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Appendix A: Analysis methods 



used to construct the smearing matris0- These three 
sources are accounted for in the covariance following the 
standard uncertainty propagation. The first two terms 
are the diagonal covariances of tT[fj7 (the raw yields), 
and <Ju^^^'^ (the background yields). Together, these give 
the uncertainty of background-subtracted yields, i.e., of 
(jdata = o-raw _ ^^1%''^^^ ^ rpj^g ^^liid term accounts for 
the propagation of the smearing matrix covariance C's 
through the full unfolding and fitting procedure. Cg 
is calculated from the statistical uncertainty of the mi- 
gration matrix and the Born-level simulated yields used 
to calculate the smearing matrix. This uncertainty con- 
tributes an additive term in C of the form 



^unf ^ 



unf 

uu 



(Al) 



where cr™^ is the unfolded (Born-level) yield vector, 
which is calculated by correcting the measured yields for 
smearing and background. This third contribution to the 
covariance is small, as the Monte Carlo productions con- 
tain approximately 20 times as many events as the data 
productions. 

The defined in Eq. llOl was minimized with respect to 
the vector of parameters /3 by means of linear regression, 
producing the parameters 



T/^-l data 



/3={X' S' C-'SX)-'X' S' C 
along with their covariance 

Cp ={X^S'^C'^SX)-^. 



(A2) 



(A3) 



In this section some technical details particular to this 
analysis are discussed. Section lA II includes the details 
of the fully differential unfolding and fitting procedure. 
Section 1X2] describes the extraction of a AD model of the 
measured cosine modulations from the data. 



Here, a generalized procedure was used that unfolds sep- 
arately each data set with its own smearing matrix, and 
then simultaneously fits the unfolded yields from the var- 
ious year^. For this purpose, a super matrix form of 
Eq. [To] was defined. 



1. Five-dimensional unfolding and fitting 



As described in section BlI B[ the measured yields are 
simultaneously unfolded (i.e., corrected for acceptance, 
smearing and QED radiative effects) and fit by minimiz- 
ing the in Eq. [TO] In Equation [TOl cr'*^*^ is a vector 
and S and C are square matrices, all of dimension of the 
number of bins (5*5*6*6* 12 = 10800, see Table UU. 
The results vector /3 contains the three fit parameters 
(A, B, and C, see Eq. [9]) for each (x, y, z, Ph±) bin and 
therefore is of dimension (5 *5*6* 6) *3 = 2700. The 
product Xp gives the fit function of Eq. [9] in each of the 
10800 bins and so X is 10800 by 2700. Each row contains 
elements equal to 1, cos(0), and cos2((/)) for that bin; the 
remaining elements are 0. The result is that X is block 
diagonal, with blocks of dimension 12x3. 

The covariance C includes the sum of three sources 
of statistical uncertainties: the precision of the measured 
yields, the precision of the Monte Carlo used for the back- 
ground subtraction and the precision of the Monte Carlo 



= {rrl^ - S^Xl3fC^\cTl^ - S«X/3), (A4) 

where cr^^j is a super vector that includes the aJu'^ ■ for 
all data sets d = 1, .., n 



cr,„, — 



data.l 
O'uu 



data.n 
CTuu ' 



(A5) 



The super matrices and Cs include respectively the 
smearing matrices and their covariances for various data 



^ It can be shown mathematically that if the results are first un- 
folded (including only the first two contributions), and then fit 
(including the additional uncertainty due to S), an identical term 
will appear in the calculation of the covariance of the fit param- 
eters. 

* Mathematically this is equivalent to taking the weighted average 
of the fit parameters from various data sets before calculating the 
moments (which correspond to the ratio of the fit parameters). 
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sets: 




Cf^ 
.. 




(A6) 



Since Cs is a block-diagonal matrix, equation I A4I gives 
X^ir^SjC^'atT''^ {AT) 

and covariance 

'^P = SjC^'S^j x| . (A8) 

2. 4D model of the extracted moments 

Several of the systematic tests described in section IIVI 
require a Monte Carlo production that includes az- 
imuthal modulations consistent with those found in the 
data. To facilitate this, the fully differential final results 
are fit to a four-dimensional function, which is then used 
to alter the underlying distribution in an azimuthally in- 
dependent PythiaG Monte Carlo production. The fit 
function used has 38 parameters, 19 for each modulation 
(cos 20, cos (/)) , and is of the form: 

/ =Ai + A2X + A^y + Aiz + A^Ph^ 
+ Aqx^ + Atz^ + AgPhi^ + A^xz 
+ AioxPhi_ + AiizPhi_ + Ai2yPh± 
+ A^yx + Auyz + Ai^y'^ 
+ Aidx"^ + Anz^ + Ai8P,,_L^ + Ai9y\ 



(A9) 



Several other functional forms were also tested and gave 
compatible results. This procedure was used to extract 
one model separately for hydrogen and deuterium tar- 
gets, and for each particle type and charge: pious, kaons, 
and unidentified hadrons. 



Appendix B: Hadron Identification 

The Hermes dual-radiator ring-imaging Cherenkov 
(RICH) detector is described in detail in Ref. (sO]. In 
that article the indirect ray tracing (IRT) particle iden- 
tification algorithm is presented. In addition, an alter- 
native method, the direct ray tracing (D RT) algorithm 
was developed, which is described in Ref. [83]. The IRT 
algorithm calculates an expected Cherenkov angle and 
compares this to the observed photons. The DRT al- 
gorithm performs a simulation, generating an expected 



photon pattern for a given particle hypothesis, which is 
then compared to the observed pattern. Because DRT 
performs a full simulation it is more accurate than IRT, 
at the cost of increased computing time. 

Here, a new method, EVT, is presented. This event- 
level algorithm can more effectively identify tracks in 
semi-inclusive events where rings from several tracks may 
overlap. 



1. The EVT algorithm 

The DRT algorithm generates a simulated photomul- 
tiplier tube (PMT) hit pattern for each radiator, based 
on known track kinematics and a particle-type hypothe- 
sis. Hit patterns are generated for each particle type hy- 
pothesis (pion, kaon, proton) and the likelihood of each 
hypothesis is calculated by comparing the simulated hit 
pattern to the hit distribution observed. Due to comput- 
ing constraints, DRT was initially only used on a subset 
of particle tracks, and the event-level sum over tracks 
shown in the equations of section 3.1 of Ref. [ssj was not 
computed. Advances in CPU power in recent years made 
it possible for the DRT method to be run on all tracks in 
all events, but the event-level consideration of all tracks 
in an event was not implemented in the software. 

The EVT method implements the full event-level iden- 
tification algorithm described in Ref. (ssj. Simulated hit 
patterns for each track are combined for all permuta- 
tions of particle hypotheses to form a set of predicted 
hit patterns for the event. The likelihood {L^) of an 
event hypothesis H is given by the observed hit pattern 

(Cpmt(»)) and the probability of hit (^'pmt(*))' 
L^ = ^log [pAitWC^pmt(«) 

i 

+ i>?MTW(l-CpMT(z))], (Bl) 

where the probability of no hit is simply ^'pmt(*) ~ ^ ^ 
-fpMT(«)- The observed hit pattern CpMT(i) is 1 if PMT 
i is hit, and otherwise. An event hypothesis H is a 
set of particle type hypotheses, one for each track in the 
event. Given T tracks, each with h (=3, pious, kaons, 
proton) possible hypotheses, there arc a total of /i^ event 
hypotheses H. The hypothesis for track t given the event 
hypothesis H is Hf. The simulated hit pattern of track 
t, of particle type Ht, produced from radiator r is given 
by Af('f^*'*''^)(i). The probability Ppuri^) of a hit given 
the event hypothesis H, is 




E 



1 



JHt,t,r} 



B{^) 



(B2) 



where the sum is over the hits from both radiators 
(r) and all the tracks (t) in the event. An unphys- 
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ically high number of hits is simulated for each ra- 
diator {J2i^^^'"*''^H^) = 360) to construct a smooth 
distribution of the expected hits. The simulated hit 
pattern(Af(^*'*''')(i)) is then normalized ( ^ '^(JJt.C)(-;) ) 
to the number of expected PMT hits for the given particle 
type, track kinematics, and radiator. The total number 
of expected hits, n^^*'*^''), is typically 0-10 hits. The B{i) 
term is included to take into account physical and exper- 
imental backgrounds; see section IB^ for more details. 

After the likelihoods are computed, the most likely 
is chosen and the particle type of each track in the event 
is given by Ht. 



2. Backgrounds 

The background term B{i) was investigated by count- 
ing the average number of hits in each PMT in the 
absence of tracks in that detector half. Average back- 
ground values of 0.005 were observed in the data, with 
a structure that shows the highest values nearest to the 
beam, indicating that the origin was not only electronic 
noise. Event displays showed that these "background" 
events that have no associated tracks in fact show a 
ring structure in the RICH. Extracting the same quan- 
tity from Monte Carlo simulation showed similar results, 
confirming that these "background" rings are from phys- 



ical particles not tracked by the spectrometer. They 
are caused by high-energy photons from tt^ decays and 
bremsstrahlung of the scattered beam lepton, producing 
e+e~ pairs in front of the RICH. 

While the B{i) term can only provide an average treat- 
ment of such untracked rings, it is the simplest way to 
take this physical background into account using the ex- 
isting algorithm. Monte Carlo tests using B{i) deter- 
mined from data showed that it led to more efficient par- 
ticle identification with less contamination compared to 
the constant value of 0.0001 (as was used in the past). 

Extracting B{i) from various subsets of the data 
showed that B{i) is relatively insensitive to the type 
of target gas. However, it is affected by the gas den- 
sity, showing higher values for targets with a higher den- 
sity. When applied to simulated data, a larger B{i) 
favors identification (both correctly and incorrectly) as 
a (anti)proton while a smaller B{i) favors identification 
(both correctly and incorrectly) as a kaon. In each case 
there is a trade off between efficiency and contamination. 
Since the overall flux of pions is largest, it is relatively 
unaffected by the change while the very small antiproton 
flux shows relatively large changes when using a different 
B{i). However, the use of different background files had 
a negligible affect on the cosine modulations extracted 
here. B(i) was extracted from the unpolarized data for 
each data year and used for the RICH algorithm applied 
to all of the data from that year. 
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